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Abstract 

The DIS diffractive cross section, da'^I^^xN /^^x , has been measured in tlie mass range 
Mx < 15 GeV for j*p cm. energies 60 < < 200 GeV and photon virtuahties = 7 
to 140 GeV^. For fixed and Mx, the diffractive cross section rises rapidly with W, 
da^!l^j^j^{Mx,W,Q'^)/dMx oc W^'"' with a'^'ff = 0.507 ± OmA{stat)tofi{syst) cor- 
responding to a t-averaged pomeron trajectory of = 1.127 ± 0.009(sfat)to oi2('5y'St) 
which is larger than cF^ observed in hadron-hadron scattering. The W dependence of the 
diffractive cross section is found to be the same as that of the total cross section for scat- 
tering of virtual photons on protons. The data are consistent with the assumption that 
the diffractive structure function F^^^^ factorizes according to XjpF^^^\xjp, /SjQ'^) = 
{xo/xjp)'^F^^'^\P,Q'^). They are also consistent with QCD based models which incor- 



porate factorization breaking. The rise of x^F^^^^ with decreasing Xjp and the weak 
dependence of F^^'^^ on suggest a substantial contribution from partonic interactions. 
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1 Introduction 



Diffraction has been studied extensively in hadron-hadron scattering at small momentum trans- 
fers [||. An elegant parametrization of the data has been provided by the Regge formalism 
through the introduction of a pomeron trajectory - The hypothesis that diffraction may 
have a partonic component has been substantiated by the observation of a high transverse 
energy jet production in diffractive pp scattering However, in hadron-hadron scattering 
both collision partners are extended objects which makes the extraction of the underlying par- 
tonic process(es) difficult. In deep-inelastic electron-proton scattering (DIS), on the other hand, 
the virtual photon has a pointlike coupling to quarks. The ep collider HERA offers a unique 
opportunity to study the partonic structure of diffraction since it gives access to the regime 
of large photon virtualities (Q^ = 10 - 1000 GeV^) and large energy transfers between the 
virtual photon and the target proton in its rest system, u = j {2mpx) = 2-20 TeV, where x 
is the Bjorken scaling variable and rup is the proton mass. 

The diffractive dissociation of the virtual photon, first recognized by the presence of a class of 
events with a large rapidity gap p, |[ , has opened a window for a systematic study of diffraction 
in reactions initiated by a hard probe [IU| - |]T3|. 



In this paper, we present a measurement of the diffractive cross section for 

7*p XN (1) 

?D(3) 



and of the diffractive structure function F2 [0. Here, X and N are the particle systems 



produced by dissociation of the virtual photon and the proton. The measurements show that 
diffraction constitutes a substantial fraction of the total cross section. The latter is directly 
related to the proton structure function F2{x,Q'^). The principal signatures for a partonic 
behaviour in DIS at small x were found to be a logarithmic dependence of F2 on associated 
with a rapid rise as x decreases [|^, |1^, |19]. The comparison of F-^^^^ with F2{x,Q'^) allows a 



2 



direct comparison of the QCD evolution of the two processes with respect to x and Q 

In QCD, diffraction is characterized by the exchange of a colourless object, e.g. a colour 
singlet two-gluon system, between the incoming virtual photon and proton. The exchange of 
a colourless system suppresses QCD radiation, and therefore the production of hadrons, in 
comparison with nondiffr active scattering. In the diffractive events studied in this analysis 
most of the hadronic energy is carried away by a low mass nucleonic system which escapes 
detection. This property is used to identify the diffractive contribution. The diffractive cross 
section is determined using the Mx method developed previously to separate the diffractive 
and nondiffractive contributions 



The present measurement is based on a fivefold larger data sample in a wider range in Q"^ (7 - 
140 GeV^) compared to our previous studies [jlT], |T3|. The squared momentum transfer t from 



the virtual photon to the incoming proton is not measured, so the diffractive contribution was 
integrated over this variable. The system A^ is either a proton or a nucleonic system with mass 
Mjy < 5.5 GeV. The new results supersede those presented in [|13] which were affected by a 
technical error simulating the QED radiative corrections in the Monte Carlo generation used 
for unfolding. This led to a steeper energy dependence and a higher intercept of the pomeron 
trajectory by about one unit of the quoted error. 
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2 Kinematics 



The kinematic quantities used for the description of inclusive DIS, e{k) + p{P) e(fc') + 
anything, are = -q^ = -{k - k'f, x = Q^/{2P ■ q), y = (P ■ q)/{P ■ k) and = 
— x)/x + rrip ^ Q^jx for x <^ 1. Here k,k' are the four- moment a of the initial and 
final state positrons; P is the four-momentum of the intial state proton and y is the fractional 
energy transfer to the proton in its rest frame. For the range of and W considered in this 
paper W"^ ~ ys, where s = iEgEp is the square of the ep cm. energy, = 300 GeV. The 
scaling variables used to describe DIS diffraction are given by x^ = [{P — N) ■ q]/{P ■ q) ~ 
(Ml- + Q^)/{W^ + Q2) and /3 = QV[2(P - N) ■ q] = x/x^ ^ QViM^ + Q^) where N is 
the four-momentum of the outgoing nucleonic system and Mx is the mass of the system into 
which the virtual photon dissociated. In models where diffraction is described by the t-channel 
exchange of a system, for example the pomeron, x^ is the momentum fraction of the proton 
carried by this system and (3 is the momentum fraction of the struck quark within this system. 



3 Experimental conditions 

The experimental conditions in 1994 for HERA and the ZEUS detector were described in our 
previous paper dealing with the F2 measurement HERA operated with 153 colliding 

bunches of 27.5 GeV positrons and 820 GeV protons. Additional unpaired positron (15) and 
proton (17) bunches circulated, which were used to determine beam related background. The 
data of this analysis corresponds to a luminosity of 2.61±0.04 pb~^. 



The ZEUS apparatus is described in detail elsewhere [gO|. Of particular importance for this 
analysis were the uranium-scintillator calorimeter (GAL) the central tracking detector 
(GTD) the small angle rear tracking detector (SRTD) |2^, the proton remnant tagger 
(PRT) |4[ and the luminosity monitor (LUMI) 



The GAL provides an angular coverage of 99.7% of An and is divided into three parts (FGAL, 
BGAL, RGAL), covering the forward (proton direction), central and rear regions with pseudora- 
pidity ranges of 4.3 < r/ < 1.1, 1.1 < r/ < —0.75 and —0.75 < t] < —3.8, respectively^. Each part 
consists of towers which are longitudinally subdivided into electromagnetic (EMG) and hadronic 
(HAG) readout cells. In test beam measurements, energy resolutions of as/E = 0.18 /\/E for 
electrons and aE = 0.35/ ^/E for hadrons were obtained {E in GeV). 

The GTD is a cylindrical drift chamber situated inside a superconducting solenoid which pro- 
vides a 1.43 T field. It consists of 72 cylindrical layers covering the polar regions 15° < 6 < 164° 
and the radial range 18.2 - 79.4 cm. The transverse momentum resolution for tracks traversing 



all GTD layers is cr{pt)/pt ~ y (0.005pt)^ + (0.016)^, with pt in GeV. The vertex position of a 
typical multiparticle event is determined from the tracks to an accuracy of typically ±1 mm in 
the X, Y plane and ±4 mm in Z. 

The PRT is used to tag diffractive events where the proton dissociated. The PRT consists 
of two layers of scintillation counters installed perpendicular to the beam at Z = 5.15 m. 



^The ZEUS coordinates form a right-handed system with positive-Z in the proton beam direction and a 
horizontal X-axis pointing towards the centre of HERA. The nominal interaction point is a,t X — Y = Z = 0. 
The pseudorapidity rj is defined as — ln(tan|), where the polar angle 9 is taken with respect to the proton beam 
direction from the nominal point. 



i.e. downstream of FCAL and beam collimator C4. The two layers are separated by a 2 
mm thick lead absorber. Each layer is split into two halves with two counters each which are 
independently read out by photomultipliers. The counters have an active area of dimensions 
30 X 26 cm^ with a hole of 6.0 x 4.5 cm^ at the centre to accomodate the HERA beam pipe. 
The pseudorapidity covered by the PRT is 4 ^ rj ^ 6. The data with useful PRT information 
correspond to an integrated luminosity of 0.7 pb^^. 



4 Reconstruction of kinematic variables 

The kinematic variables x,Q'^,W and Mx were determined from calorimeter and tracking 
information. The calorimeter cells were required to have energy deposits above 60 MeV in the 
EMC section and 110 MeV in the HAC section unless these energy deposits were isolated in 
which case thresholds of 120 MeV (160 MeV) in the EMC (HAC) sections were used. An energy- 
momentum vector {Ej,pxj,PYj,Pzj) with Ej = Vxj+Vyj'^Pzj assigned to every calorimeter 
cell j with energy deposition Ej. The cell angles were calculated from the geometric centre of the 
cell and the vertex position of the event. The algorithm used to identify the scattered positron 



was based on a neural network p6[ which included information from the CAL. The systematic 
uncertainty in the energy determination of the scattered positron is 2% at 10 GeV decreasing 
linearly to 1% at 27.5 GeV for the region considered in this analysis. The momenta of the 
particles of the system X were reconstructed from clusters found in the calorimeter and from 



tracks in the CTD |T5|, The inclusion of tracking information improves the Mx resolution 
and reduces the sensitivity to the losses due to inactive material in front of the calorimeter. A 
systematic uncertainty of 3% is assigned to the hadronic energy measurement. The selected 
clusters and tracks are called energy flow objects (EFO's). 

The kinematic variables x, and W were determined with the double angle (DA) method 
in which only the angles of the scattered positron {6'^ and the produced hadronic system (7^;^) 
are used. The angles were determined from the EFO's. In the DA method, in order that 7^^ 
be well measured, it is necessary to require a minimum of hadronic activity in the calorimeter 
away from the forward direction. A suitable quantity for this purpose is the hadronic estimator 
of the variable y ||29|, defined by yjB = T.j {Ej — pzj) /2i?e, where the sum runs over all EFO's 



not assigned to the scattered positron. 

We study events of the type ep —>■ e + X + rest, where X denotes the hadronic system observed 
in the central detector (CAL and CTD) and rest the particle system escaping detection through 
the beam holes. The mass Mx of the system X was determined by summing over all EFO's 
not assigned to the scattered positron: 

iMrn' = iEEjr~ij:px,r-iEpy^r-ij:pz,Y ■ (2) 



5 Trigger and event selection 

The event selection at the trigger level was identical to that used for the F2 analysis of the same 
data H^. The off-line cuts were also similar to those applied previously. The energy of the 
scattered positron had to satisfy E'^ > 10 GeV to ensure reliable positron finding and to suppress 
the photoproduction background. The variable y, calculated from the scattered positron, was 



required to satisfy Ue < 0.95 to suppress events with spurious low energy positrons. The 
impact point of the positron on the face of the RCAL had to he outside a square of side 26 
cm centered on the beam axis (boxcut) to ensure full containment of the positron shower. 
The requirement yjB > 0.02 ensured a good measurement of the angle jh and hence of x. 
By requiring 40 < 5 < 70 GeV, where 6 = — pzj) summed over all EFO's, including 

those belonging to the scattered positron, both the background from photoproduction and the 
radiative corrections were reduced. The primary event vertex was determined from the tracks. 
If no tracking information was present (9.2% of the events) the vertex position was set to the 
nominal interaction point. 

After the selection cuts and removal of events from QED Compton scattering and cosmic rays, 
the DIS sample contained 304k events. For the present analysis, 157k events with 60 < W < 200 
GeV, 7 < < 140 GeV^ were used. 



6 Monte Carlo simulation 



Monte Carlo simulations were used for testing the validity of the subtraction of the nondiffrac- 
tive contribution, for understanding the contribution from double dissociation (7*^ — > XN), 
for unfolding the produced event distributions from the measured ones, for determining the 
acceptance and for estimating the systematic uncertainties. The detector simulation is based 
on the GEANT program 130] and incorporates our present understanding of the detector and 



the trigger and test beam results. 

Hadronic final states from diffractive DIS interactions where the proton does not dissociate, 
ep — > eXp, were modelled with RAPGAP [31, 27| modified to include low-mass vector meson 
production. RAPGAP is based on a factorizable model in which the incoming proton emits 
a pomeron. The interaction of the virtual photon with this pomeron is described by an effective 
structure function F^{j3,Q'^) that is independent of the process of emission of the pomeron 
and where the partons of the pomeron take part in the hard scattering. The parton densities 
of the pomeron are evolved from a starting scale Ql = A GeV^ using the next-to-leading order 
DGLAP equations [0. The free parameters are adjusted to reproduce the results on the 
diffractive structure function Fg^^^"* measured by HI |T^. The momentum sum rule was not 
imposed. RAPGAP was used with the parton showering scheme of ARIADNE 4.08 ||33|, which 
is based on the color dipole model and includes the first order QCD matrix elements, and the 
Lund fragmentation scheme @] as implemented in JETSET 7.4 The region of low masses 
(below 1.1 GeV) is tuned to the measured ratio of p : resonance production [^, ^ ^ and 
allowing for a contribution from uj production. 



The contribution from the diffractive process where the proton dissociates, 7*^ XN' 



dissoc 



was 



simulated using EPSOFT [^. Assuming factorisation and a Triple Regge formalism |^ 
for modelling 7*p — >■ Xp and 7*p — > x N'^^^''^°'^ , the measured cross sections for elastic and single 



diffractive dissociation in pp (and pp) scattering, 
relate 7*p — > XN"^'"""" to YP ^P- 



pp 



pp and pp — > pN 



dissoc 



were used to 



For testing the procedure used to separate the diffractive and nondiffractive contributions, 
events from standard nondiffractive DIS processes were produced within the framework pro- 
vided by DJANGO 6.0 ^Ij. First order electroweak corrections were generated with HERA- 
CLES 4.5 [0. The positron-proton kinematics was produced with LEPTO 6.5 |Q which was 



interfaced to ARIADNE 4.08 for the simulation of the parton shower process. For fragmenta- 
tion JETSET 7.4 was used. The parton densities of the proton were chosen to be the MRSA 
set M . 



All Monte Carlo events were passed through the standard ZEUS detector and trigger simulations 
as well as the event reconstruction package. 

7 Precision of kinematic variables and binning 

The resolutions expected for the kinematical variables were estimated from Monte Carlo (MC) 
simulation. The intervals in W were chosen with equidistant bins in InW^^, thereby providing 
approximately equal numbers of events in each W bin. Here, and in the following, masses and 
energies are given in units of GeV. The average r.m.s. resolution a{lnW^) is 0.23 in the lowest 
{W,Q^) bin decreasing to 0.09 in the highest {W,Q'^) bin. A bin width of AliaW^ = 0.4 was 
chosen. For the r.m.s. resolution is less than 1 GeV^ in the lowest bin increasing to 
3 GeV^ in the highest interval. 

The mass M^*^"*, determined from the EFO's, has to be corrected for energy losses in the 
material in front of the calorimeter and for acceptance. The correction was determined by 
comparing, for Monte Carlo (MC) generated events, the MC measured mass, M^'-^"^'^"'^ , to 
the generated mass, M^'"^'^"', of the system X. The resulting corrections to determine the 
diffractive cross section were performed in three steps. 

In the first step an overall mass correction factor was calculated from the ratios of measured 
to generated masses, f{M^^~) = m^^^~ /M^^^^'^, as a function of M^, W and Q^. The 
variation of f{M^^'^^''') with M^, W, was found to be sufficiently small (< 6%) for the Mx 
range used in this analysis, 1.5 < Mx < 15 GeV, so that it could be neglected in the first step 
of the mass correction. The average correction factor was / = 0.80. The same correction factor 
was used for masses below 1.5 GeV. The correction factor / was applied to obtain the corrected 
mass value, M^°'' = M"^^""' /J . The r.m.s resolution for Mx was a{Mx)/\fM^ ~ 60%GeV3 on 
average. All Mx distributions presented below refer to Mx '^. 

For diffractive production, a comparison of the MC-generated distributions with the MC- 
measured distributions show a depletion of events at the high mass end and an excess at 
somewhat lower mass values. This (small) distortion r = jsj^^Cmeas jj^MCgen caused mainly 
by particles lost through the forward beam hole. Since r was found to be independent of Q^, 
within errors, it was determined in bins of In M'j^ for the different W intervals averaged over Q"^ 
yielding r{\n.Mx-,W)] no smoothing was applied to r. The distortion was taken into account 
in extracting the diffractive contribution. 

In the final step, the diffractive cross section was determined by an unfolding procedure dis- 
cussed below taking into account, for each {Mx-, W, Q"^) interval, the proper mass correction as 
determined from the MC simulation. 

Results are presented for the intervals in W: 60 - 74, 74 - 90, 90 - 110, 110 - 134, 134 - 164, 
164 - 200 GeV; in Q^: 7 - 10, 10 - 20, 20 - 40, 40 - 140 GeV^, with average values of 8, 14, 
27, 60 GeV^; in Mx: < 3, 3 - 7.5, 7.5 - 15 GeV, with average Mx values of 2, 5 and 11 GeV. 



8 Characteristics of the Mx distributions 



The method used to separate the diffractive and nondiffractive contributions is based on their 
very different Mx distributions. The mass distributions for typical (W, Q'^) intervals are pre- 
sented in Fig. 0(top) in terms of Mx- The distributions shown were not corrected for acceptance. 
Two distinct groups of events are observed, one peaking at low Mx values, the other at high 
Mx values. While the position of the low mass peak is independent of W , the high mass peak 
moves to higher values as W increases. 

Most of the events in the low mass peak exhibit the large rapidity gap that is characteristic 
of diffractive production. This may be seen from the shaded areas in Fig. |l] which show the 
distributions of events with rjmax < 1-5 where rjmax is the largest pseudorapidity at which 
energy deposition in the calorimeter or a track was observed. This value corresponds to a 
pseudorapidity gap in the detector larger than ~ 2.5 units since the beam hole of the FCAL is 

at r]edge - 4. 

In Fig. |I](bottom) the mass distributions are presented in terms of In M^- In this representation, 
the Mx distributions exhibit a simple scaling behaviour. The low mass peak shows up as a 
plateau-like structure at low lnM|^, most notably at high W values. The shaded histograms 
show again the distributions of the events with a large rapidity gap, rjmax < 1-5, which account 
for most of the observed plateau Q. The high mass peak has a steep exponential fall-off, 
dM/dlnMx oc exp(61nM|^), towards smaller lnM|^ values. The peak position of the high 
mass events changes proportionally to In W^, i.e. shows scaling in hi{Mx/W^) and the slope, 6, 
of the exponential in lnM|^ is approximately independent of W and Q^. These characteristics 
are properties of events with uniform, random and uncorrelated particle production along the 
rapidity axis where particles are accepted in a limited range of rapidity. In models such as the 
Feynman gas model or one dominated by longitudinal phase space [^, the slope h represents 



the particle multiplicity per unit of rapidity. ^ The exponential in InM'j^ and the scaling in 
ln(M^/l^^) are directly connected to the exponential suppression of large rapidity gaps by 
QCD radiation. The latter populates the region between the struck quark and the coloured 
proton remnant. 

These characteristics are also properties of realistic models for particle production in deep 
inelastic scattering. ARIADNE, which gives a good description of particle production by DIS 
at HERA, also exhibits a pure exponential fall-off with lnM|^ and scaling in InM^/iy^.This is 
shown in Fig. 0(top) which presents the MC simulation of the nondiffractive contribution at the 
generator level as a function of ln(M|-) for different {W, Q^) intervals (solid histograms) where 
only particles generated with i] < 4.3 were included. The slope b is independent of W and 
Q^, which is in agreement with the experimental observation that, for DIS events, the average 



charged particle multiplicity per unit of pseudorapidity t] is independent of W and Q ^ 



The comparison of the distributions at the generator level (dotted histograms) and the detector 
level (dashed histograms) in Fig. |^(top) shows that the exponential fall-off with In M"^ is not 
affected by detector effects. 

^Note that a cut on rimax will select events produced by colourless exchange, from both diffraction and from 
reggeon exchange but does not, in general, allow the extraction of the full diffractive contribution (see also p5[|). 
■^The pseudorapidity gap At] is related to InAf^ via A77 w ln(iy^/Af|-) as discussed in detail in fl^ . 



9 Extraction of the difFractive contribution 



In diffractive events, the system X resulting from the dissociation of the virtual photon is 
almost fully contained in the detector while the outgoing proton or low mass nucleonic system 
escapes through the forward beam hole. Furthermore, diffractive dissociation prefers small Mx 
values and leads to an event distribution of the form dJ\f/dMx oc 1/(M|^)*^^"'""^ corresponding 
to dJV /dlnMx oc 1/(M|^)", approximately independent of W. At high energies and for large 
Mx, one expects n ^ 0, leading to a constant distribution in In Mj^. Such a mass dependence is 
seen in diffractive dissociation of pp and pp scattering ^ A value of ^ is also expected 



in some models |Q - for diffractive DIS as a limiting value for the fall-off of the mass 
distribution. 

We identify the diffractive contribution as the excess of events at small Mx above the expo- 
nential fall of the nondiffr active contribution in InMx. This will be referred to as the Mx 
method for the determination of the diffractive component. The exponential fall permits the 
subtraction of the nondiffractive contribution and, therefore, the extraction of the diffractive 
contribution without assuming the precise Mx dependence of the latter. The distribution is 
expected to be of the form: 

D + c exp(6 In M|.), for In Mj- <lnW^ - tiq. (3) 



d In Ml 



Here, D denotes the diffractive contribution and the second term the nondiffractive contribu- 
tion. The diffractive term D is multiplied by the distortion function r{\nMx,W) discussed 
above. The quantity (In W'^ — rjo) specifies the maximum value of In Mx up to which the ex- 
ponential behaviour of the nondiffractive part holds. A value of r/o = 3.0 was found from the 
data p. We apply Eq. |^ in a limited range of InMx to determine the parameters b and c 
of the nondiffractive contribution. The diffractive contribution is not taken from the fit result 
for D but is determined by subtracting from the observed number of events the nondiffractive 
contribution found from the fit values of b and c. 



The diffractive contribution in various models of diffraction |^, is expected to be a 
slowly varying function of lnM|^ when Mx > and to approach, in the asymptotic limit, a 
constant InMx distribution at large Mx, D = constant. The nondiffractive contribution in 
the {Mx,W,Q'^) bins, was determined in two steps. This procedure will be referred to as the 
nominal analysis. In the first step, the slope b was determined as an average over the fits to the 
data of the high W intervals 134 - 164, 164 - 200 GeV for = 7-10 and 10 - 20 GeV^ in the 
restricted Mx range. In lOQ^ < lnM|. < \nW^ - r/o. The fits yielded b = b^om = 1-72 ± 0.07. 
This value is smaller than that predicted with ARIADNE by about 10% (see below). 

In the second step, the fits were repeated for all (W, Q"^) intervals in the wide Mx range, 
InQ^ < lnM|^ < In 14^^ — ?7o, using b = bnom as fixed parameter and assuming D = constant. 
The fit results for the nondiffractive contribution are shown in Fig. |I] by the straight lines, those 
for the sum of the diffractive and nondiffractive distributions by the upper curves; the latter 
fluctuate since no smoothing was applied for the distortion correction r. The fit results for 
the sum of the diffractive and nondiffractive contributions give a good description of the mea- 
sured event distributions. Subtraction of the nondiffractive contribution yielded the number of 
diffractive events. From these, the diffractive cross section was obtained by unfolding. Fits were 
also performed with different forms for D. These were used for estimating the corresponding 
systematic uncertainties (see below). 



The fit procedure for separating tlie diffractive from the nondiffr active contribution was tested 
with MC simulated event samples from diffractive (using RAPGAP) and nondiffractive scat- 
tering (using ARIADNE). It was first checked that the slope of the exponential fall-off for the 
nondiffractive contribution is not affected by detector effects. The lnM|^ distributions pre- 
dicted for the nondiffractive contribution are shown in Fig. 0(top) at the generator (dashed 
histograms) and detector (solid histograms) levels for the W intervals and values. The 
straight lines show the fits to the MC measured distributions. The value of the slope Jj'^'^aen 
at the generator level (average value Jj^^'^aen _ -^ g^ _|_ q_q3^ agrees well with the value h'^^"^^"'^ 
at the detector level (average value l)^'^^^"-^ = i.gg ± 0.03) for each {W,Q'^) interval. Fits 
performed to the sum of the diffractive and nondiffractive contributions in the restricted Mx 
range are shown in Fig. ||(bottom). They resulted in an average slope, = 2.04±0.09, which 
agrees well with the 6-values found from the fits to the nondiffractive sample alone. In this 
figure the distributions of In M'^ are displayed for the sum of the diffractive and nondiffractive 
contributions as predicted at the detector level (points with error bars) and for the nondiffrac- 
tive contribution alone (dashed histograms). In the next step, the sum of the diffractive and 
nondiffractive contributions were fitted in the wide Mx range using b = h^^. The straight lines 
show the nondiffractive contribution as obtained from the fits; they give a good description of 
the input distributions for the nondiffractive contribution shown as the dashed histograms. The 
upper curves show the fit results for the sum of the diffractive and nondiffractive contributions. 
A good description of the input distributions is observed. Using the fits with h = 6^^, the non- 
diffractive contribution was subtracted which yielded, for every {W, Q^) interval, the number 
of MC measured diffractive events. From this result, the number of MC produced events was 
determined by unfolding. Comparison with the generated (i.e. input) numbers of diffractive 
events showed very good agreement for all (Mx, W, Q"^) intervals. 

Before applying the fitting and unfolding procedure, the number of events observed in the 
{Mx,W,Q^) intervals were corrected for the contribution from positron (proton) beam gas 
scattering. Averaged over all events, the contamination from positron (proton)-gas scattering 
amounted to 2.1% (1.0%). For the nominal analysis, which did not require an event vertex, 
beam-gas contributions were small for Mx > 3 GeV and negligible for above 20 GeV^. 
Sizeable contributions were observed for Mx < 3 GeV when W > 90 GeV and < 20 GeV^ 
where they amounted to 14% on average for positron-gas scattering and < 10% for proton-gas 
scattering. To estimate the systematic uncertainties the analysis was also performed requiring 
an event vertex (see below). In this case the background from beam-gas scattering was negligible 
everywhere. 

The number of diffractive events produced, MprJd^ was obtained from the number of events 
determined from the fits, A/'mfL' by an inverse matrix unfolding procedure which took into 
account detector effects such as bin-to-bin migration, trigger biases and event selection cuts. It 
was checked that the diffractive model (RAPGAP) describes the energy flow as a function of 
7] for rjmax < 1-5 for all (Mx, W, Q"^) bins as well as the observed Mx distributions in the 
region Mx < 3 GeV. 



9.1 Contribution from nucleon dissociation 

An estimate of the contribution from diffractive dissociation of the proton, 7*p ^f^dissoc^ 
(double dissociation) to the diffractive cross section presented below was obtained by comparing 
the contributions for 7*p — > Xp (single dissociation) with an identified proton in the leading 
proton spectrometer LPS [|l^] and for 7*p — XA^ determined in this analysis. This led to the 



fractional contribution from double dissociation (XA^ — Xp)/XN = XN'^'"^^"'^ / XN = (31 ± 
15)%. 

The procedure for extraction of the diffractive contribution was tested at the MC level also 
for the case where an additional contribution from diffractive dissociation of the proton was 
present. The events for 'j*p — >■ XN'^^^^°'^ were simulated using EPSOFT. In Fig. |^ the lnM|- 
distributions are presented at the detector level for the diffractive contributions from 7*p — >■ 
^j\jdissoc (^dQif^Qd histograms) and for the sum of 7*p — > Xp + XN'^'^^^"^ (solid histograms). Also 
given are the nondiffractive contributions (dashed histograms) and the sum of the diffractive 
and nondiffractive contributions (points with error bars). Fits to the sum of the diffractive 
and nondiffractive contributions in the restricted Mx range yielded for the slope a value of 
l^MCipN) ^ ^ ± 0.08 in good agreement with the value b^^^"^""' = 1.96 ± 0.03 obtained before 
from a fit to the nondiffractive sample alone. The solid straight lines show the nondiffractive 
contribution as obtained from the fits; as before, they give a good description of the input 
distributions for the nondiffractive contribution shown as the dashed histograms and are in close 
agreement with the fits performed to the nondiffractive contributions alone (dashed straight 
lines). Using the fits and subtracting the nondiffractive contribution yielded for every (VF,Q^) 
interval the number of MC measured diffractive events; these were found to be in very good 
agreement with the number of events in the diffractive sample. This is a confirmation that 
the lnM|^ method for extracting the diffractive contribution gives reliable results also in the 
presence of diffractive contributions where the proton dissociated. 

Only a limited mass (Mat) range of ]\f<^^^^°'^ contributes to the diffractive cross section presented 
below. This was studied with EPSOFT. The distribution of the generated mass peaks 
at low values, M^r ^ 2 GeV. Due to the dominance of small Mjy values the secondaries from 
j\j-dissoc g^j^g strongly coUimated around the direction of the incoming proton. Analysis of the 
angular distribution of the secondaries as a function of Mat showed that for < 2 GeV 
basically no energy is deposited in the calorimeter, while for events with Mn > 6 GeV there 
are almost always secondaries which deposit energy in the calorimeter. Furthermore, events of 
the type 7*p — > XN'^^'^^"'^, where decay particles from X'^'^'^^°'^ deposit energy in the calorimeter, 
have in general a reconstructed mass which is substantially larger than the mass of X and, 
therefore, these events contribute little to the event sample selected for diffractive production 
of 7*p XN. 

As a result, for each Mx bin there is an Mx value, called M'^'^^^*, such that the number of 
events with Mx < M'^'^^^\ which migrate outside the bin is equal to those with Mx > M'^'^^^'^, 
which migrate into this bin. For the (Mx, W, Q'^) bins studied M'^^^^* was found to be 5.5 GeV 
to within ±1.5 GeV. The spread in the value of M^^'^^'^ introduces an uncertainty in the cross 
section measurements which is below the statistical uncertainty. Therefore, the diffractive cross 
sections are quoted below as cross sections for 7*p — > XN with Mx < 5.5 GeV. 

In the nominal analysis, the unfolding of the diffractive contribution was performed consider- 
ing only dissociation of the virtual photon as simulated with RAPGAP. In order to test the 
sensitivity of the unfolding to a contribution from double dissociation, the nominal unfolding 
procedure was also applied to a MC sample consisting of a mixture of single and double disso- 
ciation in the proportion given above for the measured region. The resulting cross section gave 
a good match to the cross section with Mx < 5.5 GeV as input. 

An attempt was made to estimate quantitatively the contribution from proton dissociation by 
studying those events where a signal was recorded in the PRT. The PRT registers particles 



emitted at large pseudorapidities, 77 = 4 — 6. Such particles may result from diffractive produc- 
tion where the proton dissociated or from diffractive dissociation of the virtual photon into a 
system with a large Mx-, and from the proton remnant in nondiffr active production. The points 
with error bars in Fig. ^ show for the subset of the data with useful PRT information the In Mx 
distributions for all events. The straight lines show the nondiffr active contribution as obtained 
from the fits performed in the wide Mx range. The sum of the diffractive and nondiffr active 
contributions as obtained from the fits (upper curves) are in good agreement with the data. 
The solid histograms in Fig. ^ show the distributions (called PRT-tagged distributions) for 
events with a PRT signal. Most of the events in the high lnM|^ peak have a PRT signal as 
expected for nondiffractive scattering. The exponential fall-off of the PRT-tagged distributions 
agrees well with the predictions obtained from the fits for the nondiffractive contributions as 
shown by the straight lines. This gives additional experimental proof for the reliability of the 
Mx method. 

A fraction of the events in the low mass region also has a PRT signal. MC simulation showed 
that for Mx < 3 GeV almost all events with a PRT signal are diffractive events where the 
proton dissociated. From the number of PRT tagged events in this region, using the efficiencies 
determined for the PRT counters and the prediction from EPSOFT for the PRT acceptance, 
events from diffractive dissociation via 7*p — > XN'^^^^°'^ contribute an estimated fraction of 
j^jydissoc^j^jY _ 27 ±2{stat)% (plus an unknown systematic uncertainty which depends on the 
mass spectrum and decay properties of ]\f'^^-'^^°<^^ to the diffractive cross sections presented below. 
We prefer the estimate obtained by comparison with the LPS analysis over that obtained here 
with the PRT since the former does not rely on an assumption for the properties of _/V''*^*°'=. 



10 Evaluation of the cross sections 

For the final analysis, only bins where the fraction of nondiffractive background was less than 
50% and the purity was above 30% were kept. Purity is defined as the ratio of the number of 
events generated in the bin and observed in the same bin divided by the total number of events 
observed in the bin. The average purity in the [Mx, W, Q"^) bins was 52 %. The nondiffractive 
background fraction was small for low Mx or large W values: for Mx < 3 GeV the background 
fraction was typically below 2% except for = 27 (60) GeV^ where in the lowest W bin for 
which cross sections are presented it amounted to 10% (17%). For 3 < Mx < 7.5 GeV the 
background fraction was typically below 5% except for the next-to-lowest (lowest) W bin where 
it reached 18% (33%). For 7.5 < Mx < 15 GeV the numbers were 22 - 37% in the lowest bin 
and typically 15% elsewhere. 

The average differential cross section for ep scattering, in a given (Mx, W, Q"^) bin, was obtained 
by dividing the number of unfolded events, J^prJ^^ by the luminosity, the bin widths and the 
QED radiative correction factor. The lower limit of Mx was taken to be 2m.,^, where m^r is the 
pion mass. 

The cross section for the process ep eXN can be expressed in terms of the transverse (T) 
and longitudinal (L) cross sections, a^^^ and cr^-^^, for 7*p XN as [p5|]: 



da^l^^^xAMx, W, Q') ^ dj^ff + aff) ^ daf^UxN^Mx. Q') 

dMx dMx ^ a {1 - + 1 dMxdlnW^dQ^ ' 



Here, a term [1— di//^ dt// ] multiplying [a^^^ +af^^] has been neglected. Since y ~ W'^/ s, 

this term can be substantially different from unity only at high W values. The effect is less 
than 5% ii W < 158 GeV corresponding to y ^ W^s < 0.28, or if < 200 GeV and 
cr'^^-^ < 0.5a^^^ . In the extreme case that a"}^^^ 3> ctt^^ , the term will increase [a^^^ + cr'^^^] by 
at most 11% for the highest W bin (164 - 200 GeV). 

The differential cross section da^Ip^xN/^Mx was determined from Eq.(||) in the different 
{Mx,W,Q'^) intervals using the fit results obtained assuming a constant D and transported 
to convenient values of Mx and using the shape of the parametrization of (see below). 



Systematic errors: The systematic uncertainties for the cross section were estimated in a 
similar way as in the previous work ||13[ by varying the cuts and algorithms used to select 



the events at the data and MC levels and repeating the full analysis for every variant. The 
uncertainties in selecting the DIS events arising from the identification of the scattered positron 
were estimated by varying the box cut, the minimum energy and the minimum probability 
required for the scattered positron. The uncertainties in the reconstruction of the hadronic 
system were estimated by varying the yjB cut between 0.01 and 0.03. In each case the deviations 
from the nominal result were well below 10% and typically 1 - 3%. The uncertainty due to the 
beam-gas background was determined by removing events without a reconstructed event vertex. 
This resulted in small differences of up to 13%. In order to test for remaining background from 
photoproduction and the sensitivity to radiative effects the cut on J2j{Ej — pzj) was varied. 
This resulted in small changes except for three (Mx, W, Q^) bins where differences of up to 14% 
were seen, commensurate with the corresponding statistical errors. 

Variation of the noise cuts for isolated calorimeter cells had a negligible effect. Allowing for ±5% 
differences in the Mx correction between data and MC simulation to account for uncertainties 
in the energy scale of the calorimeter and the material in front affected mainly the low Mx bin, 
where at higher values, differences of up to 11% were found. Note that increasing (lowering) 
the mass correction factor systematically decreases (increases) the cross section. Since the 
relative changes are basically independent of W this does not affect the W dependence of the 
cross section. The acceptance corrections and unfolding were also performed with a different 
model for diffraction For < 27 GeV^ the differences were typically a few percent reaching 
values between 11% and 18% in four bins. In the highest bin the differences were < 14% 
except for the low Mx bin where they reached up to 57%; this is mainly due to the small event 
numbers involved. 

The nominal fitting procedure for subtracting the nondiffractive background used the fixed 
slope b = hnomi see Section ^ The analysis was also performed with h determined for the high 
W bins 134 - 164, 164 - 200 GeV from the wide Mx range using for D the extended form: D = 
do{l-/3)[(3{l-(3) + ^{l-/3)% and the variable form: D = do (l-/3)[/5(l-/3) + f(l-/3)5], where 
do, di and g are parameters ||50| - |Q. Both forms consist of a quark-like (first term) and a gluon- 
like (second term) contribution. In the wide Mx range the quark-like contribution decreases 
with In Mx while the gluon-like term increases. In the fits with the variable form the gluon-like 
term dominates at large masses which leads to D = constant in the asymptotic limit. The 
resulting changes of the diffractive cross section were small and did not exceed 13%. For the 
extended form, typical differences of 10 - 25% were found; the largest difference was found for 
one bin at = 60 GeV^ and amounted to —37%. 

^In this case diffraction was modelled in RAPGAP using the parametrization developed in on the basis 
of the model of iBOl. 



The total systematic error for each bin was determined by quadratically adding the individual 
systematic uncertainties, separately for the positive and negative contributions. The total 
errors were obtained by adding the statistical and systematic errors in quadrature. The errors 
do not include an overall normalization uncertainty of 2% of which 1.5% is from the luminosity 
determination and 1.0% from the uncertainty in the trigger efficiency. 



11 Differential cross section for 7*p ^ X N 
11.1 W and dependence of da'!^fp_^xN/d'Mx 

The diffractive cross section for 7*j9 — XN, {Mjsf < 5.5 GeV), is presented in Table ^and Fig. |] 
as a function of W for various Mx and values. From comparison with published data [B^ 



38], about 20% of the diffractive cross section observed in the mass bin Mx < 3 GeV at 
7 < < 20 GeV^ results from the production of the vector mesons V = p°, via 7*p VN. 
In Fig. H Q'^ da'!^!p^xN / X is presented as a function of for different (Mx, W) values. A 
fast decrease with is observed for small Mx which is similar to the behaviour of DIS vector 



meson production [g^ - ||38|. The decrease becomes slower for the high Mx region showing 
that high mass dissociation of the virtual photon becomes increasingly more important as 
grows. For a discussion in terms of a partonic behaviour see below. 

The diffractive cross section rises rapidly with W at all values for the Mx bins up to 7.5 GeV. 
The cross section was fitted for each (Mx, Q"^) bin using the form 

da'!;!l^xN{Mx,W,Q^) 



dM 



X 



h-W^" , (5) 



where a'^^^^ and the normalization constants h were treated as free parameters. The form given 
m Eq. I with a single value of a'^^^f gives an acceptable value {x^/d.o.f. = 37/35 = 1.1), 
considering only the statistical errors) although one observes in Fig. ^ (see dashed curves) a 
tendency for the data at low Mx, high to have a steeper slope than resulting from this fit. 
The fit yielded a"^'^^ = 0.507 ± OmA{stat)toMlisyst). Here, the largest contribution to the 
positive systematic error arises when the extended form is used for D in the determination of 
the slope b. In Regge models, a'^*-^-^ is related to the trajectory of the pomeron a pit), averaged 
over t: = 1 + a'^'f^ /A. The fit value for a'^'f^ leads to = 1.127 ± 0. 009 {stat)toM^{syst). 

In order to test for a possible dependence of cF^ a fit was performed where a'^^^-^ was taken 
as a free parameter for every (Mx, Q'^) bin. This resulted also in a good description of the data 
(see solid curves in Fig. |^). The resulting values are shown in Fig. ^with the statistical and 
systematic uncertainties added in quadrature. The systematic uncertainties were estimated by 
repeating the fits independently for every source of systematic error. With the present accuracy 
of the data, no conclusion can be drawn on a possible dependence of a^. 

In deriving the diffractive cross section, the assumption was made that the term 

[1 ~ (1-1)2: 1 difs^ d^fs ] can be neglected (see Eq. ^). This holds when af^^^ <^ a^^^ ■ If the 

assumption is made that a'^l^^ = a^^^ {a^^^ <^ o't^'^) ^^e value of a'^*-^-^ increases by 0.047 
(0.096) and hence increases by 0.012 (0.024). 



HI IjT^ has given for the intercept of the pomeron trajectory a value of ajp{0) = 1.203 ± 
0.020{stat)±0.013{syst)t'^o!^§{model). Avera ging over the t-distribution f\ gives approximately 
cF^ = C(jp{0) — 0.03, a value which is consistent with the result from this analysis 0. 

Our value of (T^ = 1.127 ± 0. 009 {stat)1lQ^ll{syst) lies above the results deduced from hadron- 
hadron scattering where the intercept of the pomeron trajectory was found to be (0) = 

1.08 |] and a^'"^^ (0) = 1.096t 

aoo9 O- Averaging over t reduces these values by about 



0.02 [feading to a^'°f^ = 1.06 [§ and a^'°^^ = 1.076t[J:[}o9, respectively. The latter value is 
shown by the horizontal bands in Fig. 0. 



11.2 Comparison of the diffractive and total cross sections 

The ratio of the diffractive cross section to the total virtual-photon proton cross section, 



"tot = ' (6) 



IP 



is displayed in Fig. | as a function of W for the different Mx bins (M^ < Mx < Mb) and 
values. The total cross section was taken from our F2 measurements performed with the 
1994 data @ using a%{W,Q^) = Jj'^ . F2{x ^ ^,Q^). The data show that, for fixed 



Mx, the diffractive cross section possesses the same W dependence as the total cross section. 
A fit of rf^t^ using the form rf^l^ oc W^, considering all data with Mx < 7.5 GeV and 
treating p and the normalization constants for every (Mx,Q^) interval as free parameters, 
yielded p = 0.00 ±0.03(stat), consistent with W independence. The same conclusion is reached 
when comparing the value of the power a'^*-^-^ = 0.507 ± 0.034(stat)lo'j4g(s?/st) with the power 
a*°* = 0.55 ± 0.02 obtained for 0"*°^ in the same [W, Q"^) range. Equal powers for the diffractive 
and total cross sections is contrary to the naive expectation. Assuming (i) for DIS the validity 



of the optical theorem |59] and (ii) that the cross section for diffractive photon dissociation 



at fixed Mx has the same W dependence as the forward cross section for elastic scattering, 
7*p — *• 7*p, then da^f^^xN/^^x should be proportional to W°- if (7*°^ oc W"-^"^. Hence, taking 
the W dependence found for da'!^!p^xN / dMx the power p would have been expected to be 
p = a'^^^f /2 = 0.25 ± 0.02(stat)lo;o2(s?/st), in clear disagreement with the data. 



The rapid rise of atot with W , which is equivalent to the rapid rise of F2 as x — > 0, in QCD is 
attributed to the evolution of partonic processes. The observation of similar W dependences 
for the total and diffractive cross sections suggests, therefore, that diffraction in DIS receives 
sizeable contributions from hard processes. The same W dependence for the diffractive and 
total cross sections was predicted in (see below). 



The diffractive contribution to the total cross section for low Mx decreases rapidly as 
increases while for Mx = 11 GeV it is the same, within a factor of two. Averaged over W = 



^When parametrizing the ^-distribution by da/dt cx exp[(bo + 2a jp' Inll / x jp))t\ the difference a^(0) 
determined mainly by &o and a^'. HI quotes 60 = 4.6 GeV~^ and a^' = 0.26 GeV~^. 



^In the analysis of HI |14 the diffractive contribution was first extracted by selecting events with a large 
rapidity gap. From the Xjp dependence of the resulting diffractive structure function F^^^\xjp, P, Q'^) it was 
concluded that at large a substantial reggeon contribution was present. The latter was determined from a 
Regge fit to the data using the sum of a pomeron and a reggeon contribution. We have performed a similar fit 
to our data and found no evidence for a significant contribution from reggeon exchange. 

"^Assuming for the ^-distribution da/dt oc ea:p[(6o + 2a jp' hiW'^)t\ with = 7 GeV^^ (from vr^p clastic 
scattering, see §) and a^' = 0.25 GeV^^ Q 



60 — 200 GeV (considering only the accepted bins), the diffractive contribution (with Mjv < 5.5 
GeV) to the total cross section for Mx < 7.5 GeV is 9.2t[|i% (6.3lgJ%, 3.9l|jJ%, 1.3t[i:i%) at 
Q2 = 8(14, 27, 60) GeV^. The corresponding fractions for Mx < 15 GeV in the W range 134 
- 200 GeV are 13.2t|j:f% (9.4t[||%, 7.5t|]J%, 3.7t[ii%). As increases, the high Mx region 
becomes more important. 

The behaviour of the diffractive contribution to the total cross section observed here for 
Mx < 15 GeV is not in contradiction to our earlier finding that the fraction of DIS events with 
a large rapidity gap (LRG) is rather constant with Q"^ p: the LRG analysis considered events 
with masses Mx ^ 4 GeV, did not impose an upper limit on Mx and included the W range up 
to 270 GeV. 



12 Diffractive structure function of the proton 

The concept of a diffractive structure function introduced in p is based on the assumption 
that diffraction is mediated by the exchange of a colourless object, called a pomeron, which 
is composed of partons. The diffractive structure function of the proton can be related to the 
diffractive cross section as follows [|1^: 

1 da^!U.{Mx,W^Q-) r.,^ 
2Mx dMx "^"^ "(Q2 + iy2)2g2^2 (7j 

For W'^ ^ Q^, which holds for this analysis, Eq. |^ can be written as: 

1 da';!l^xN{Mx.W^Q') ^ A.^a ^s, . . . 

2Mx dMx ^Q\Q^ + Ml)^^ [P^^^^^ )■ W 

If ^2^'''^^ is interpreted in terms of quark densities then it specifies for a diffractive process the 
probability to find a quark carrying a momentum fraction x = jSXp of the proton momentum. 

It has been suggested P that F<^^^\xjp, l3,Q^) should factorize into a term which depends on 
the probability of finding a pomeron carrying a fraction x^ of the proton momentum and the 
pomeron structure function F^^'^^ given in terms of the pomeron's quark densities which depend 
on P and Q^: 

F,''^'\x,,P,Q') = f^{x,)F^^'\p,Q') (9) 
where fip{Xp) is generically called the pomeron flux factor. 

The quantity XjpF<^^^^ is given in Table |] and shown in Fig. ^ as a function of Xjp for different 
combinations of (3 and Q^. 

In Fig. the data from this analysis (solid points) are compared with ZEUS data obtained 
using the Leading Proton Spectrometer (LPS) [|l^ and with those of HI 0]. For ease of 
comparison the XjpF^^^^ values from this analysis were scaled to the (/3,Q^) values used in 
the HI analysis. The LPS data correspond to events of the type 7*p Xp with an identified 
proton. No correction was applied for the contribution from double dissociation which is present 
in this analysis but not in the LPS data. The correction would increase the LPS data by a 
factor of 1.45lo;23. There is consistency between this analysis and the LPS data. The HI data 
correspond to Mx < 1-6 GeV while those from this analysis are given for Mx < 5.5 GeV. No 
correction was applied. The data from HI approximately agree with those from this analysis. 
However, for fixed (3, the HI values have a tendency to rise faster with even allowing for an 
extra scaling factor. 



12.1 Xp dependence of XpF2 {Xp,/3,Q'^) 



The data from this analysis (Fig. ^ show that XjpF:^^^\xjp, (3, Q^) decreases with increasing x^, 
which reflects the rapid increase of the diffr active cross section with rising W. Assuming the flux 
factor to be of the form fip{Xjp) = (C/x^p) ■ (xo/Xp)^, taking for the arbitrary normalization 
constant C = 1 and for xq the average value of the measured x^, Xq = 0.0042, the data 
were fitted with the form XjpF^^^\xjp, (3,Q'^) = (C/Xp) ■ {xo/Xjp)'^F^^'^\(3,Q'^), which leads 
to F^^^\p,Q^) = xoF2^^^^(xo,/3,Q2). The values F^^^\f]i,Qf) at the 12 measured 
points and n were treated as fit parameters. A good fit was obtained {x^/d.o.f. = 41/40 = 1.0, 
statistical errors only) yielding n = 0.253 ± 0. 017 {stat)1:omUsyst). Note that n ~ a'^'^^ /2 = 
2{a^—l). The fit was also performed assuming n to depend logarithmically on Q^. This resulted 
in small differences that were included in the errors given for F^^'^l The fact that a good fit 
was found with a single value for n shows that the data are consistent with the assumption 
that ^2^''^'' factorizes into a fiux factor depending only on x^ and a structure function F^^"^"^ 
which depends on {f3, Q"^). Note that there is an arbitrary normalization factor for the fiux and 
therefore also for F^^'^K 

12.2 f] and dependence of XpF2^^\xp, /3, Q^) and F2^'^\p, Q^) 

The F.^^'^\f3,Q'^) values obtained from the fit described above are presented in Fig. ^ as a 
function of (3 for all values. It should be stressed that these F^^'^^ values do not depend 
on whether the F^^^^ factorizes into a pomeron fiux factor or not since a) F'^^'^\j3 , Q"^) = 
xoF.['^^\xo, f3,Q'^) and the fit was basically only used to interpolate to x^ = xq] b) a fit with 
a dependent fiux gave basically the same F^^"^^ values. The data show that the diffractive 
structure function F-^'^'^^ has a simple behaviour. For (3 < 0.6 and < 14 GeV^, F^^'^^ is 
approximately independent of (3. For j3 < 0.8 also the data from different values are rather 
similar suggesting a leading twist behaviour characterized by a slow InQ^ type rescaling. For 
P > 0.9 the data show a decrease with (3 oi Q"^ . 



The behaviour of XjpF^^^\xjp, f3,Q'^) is shown as solid points in Fig. 0. The data are 



presented for fixed (Mx, W), the variables in which the diffractive contribution was extracted. 
Given {Mx,W) and Q"^ the value of be calculated. For Mx < 7.5 GeV, XjpF2 

decreases with while for Mx = 11 GeV it is approximately constant. 



D(3) 



Strong variations, which are found e.g. for the diffractive cross section (see Fig. are just 
a reflection of kinematics: the strong variation of Q'^da'^Ip_^xN/d^x is mainly controlled 
by the kinematical factor Mx/iQ"^ + M|-) in Eq. ||. 

The approximate constancy of Fg^^^"* for (3 < 0.9 combined with the rapid rise of F^^^"* as x^ 
decreases can be interpreted as evidence for a substantial partonic component in DIS diffraction 
dissociation. 



12.3 Comparison with models 



The diffractive process in DIS has attracted considerable attention because of the possibility 
that this process can be described by perturbative QCD (pQCD). In parton models the process 



can be visualized as a fluctuation of the incoming virtual photon into a qq pair followed by the 
interaction of this pair with the incoming proton leading to a qq state plus, well separated in 
rapidity, a proton or debris from the dissociation of the proton. In it was argued that the 



dominant contribution to diffraction in DIS comes from the aligned jet configuration where q 
and q from photon dissociation have small transverse momenta relative to the direction of the 
virtual photon leading to the same energy dependence as observed for diffraction in hadron- 
hadron scattering. This contribution was expected to scale with Q^. The (3 distribution for the 
aligned jet configuration from transverse photons was predicted to be of the form 

F^^^oc/?(l-/3). (10) 

The same j3 dependence was expected in pQCD when the aligned quarks interact with the 
proton through two-gluon exchange |Q . The production of a qqg system by transverse photons 



was also found to be leading twist and was assumed to have a [3 dependence of the type fS^ 



F'^{i~(3r (11) 



with 7 = 2. A later calculation |Q found 7 = 3. In the same approach the contribution to the 
production of a qq system by longitudinal photons was found to be of higher twist and to have 
a (3 dependence of the form 



f;^^ oc (3\l - 2(3Y. (12) 

In pQCD models the Xjp dependence is expected to be driven by the x dependence of the square 
of the gluon momentum density of the proton [^, [x ■ (^(x, ;U^)]^, with x = x^ and /i is the 
probing scale. 

In the sum of the contributions from the three terms F^, F^g was evaluated in the 
perturbative region. 

We now compare the data with three partonic models (NZ) (BPR) and (BEKW) g. 
In the NZ model, diffractive dissociation is described as a fluctuation of the photon into a qq or 
qqg Fock state [5^. The interaction with the proton proceeds via the exchange of a BFKL [3C 



type pomeron, starting in lowest order from the exchange of a two-gluon system in a colour- 
singlet state. The BPR model describes the process 7*p Xp as the scattering of a colour 
dipole from the photon on a colour dipole from the proton. The model parameters were chosen 
by comparison with the HI data The predictions of the NZ and BPR models are shown 
in Fig. ^ for function of Xjp by the dashed and dotted curves, respectively. The 

NZ model provides a reasonable description of the data. The BPR model has some difficulties 
in reproducing the data for medium values of (3 and < 14 GeV^. 

We studied the individual contributions from the three terms in Eqs. [T^ following which 
identified them as the major contributors to the diffractive structure function. In they were 



calculated in the perturbative region and extended into the soft region. The Xjp dependence 
was assumed to be of the form (l/x^p)". The power n was allowed to be different for the 
transverse (n^) and the longitudinal (n^,) contributions. The normalizations of the three terms 
were determined from the data. This is called the BEKW model in the following: 

x^F2^(')(/3, x^,Q^) = CT ■ Fl + CL ■ F}^ + c, ■ (13) 



with 



fS = (14) 



Xjp 



^Xjp' Q2 I V4 ' 4^g2 







F4 = (^)-(«^)-ln(l + g).(l-/?)- (16) 
riT,L(Q') = 0.1 + 4^-ln[l + ln(^)]. (17) 

The three terms behave differently as a function of . Except for a possible dependence 
of the power n^, F'^- does not depend on as a result of the limited quark pt in the aligned 
configuration. The term is higher twist but the power l/Q^ is softened by a logarithmic 
factor; grows logarithmically with similar to the proton structure function F2 at low x. 

The coefficients CT,CL,Cg as well as the parameters nj^,n1 and xo,Ql were determined from 
experiment. In the fit the power 7 was also considered as a free parameter. Assuming 
Ql = 1 GeV^ and xq = 0.0042 and treating the other constants as free parameters a good 
fit {x^/d.o.f. = 56/47 = 1.2, statistical errors only) was obtained for the XjpF^''^\xjp, f3,Q'^) 
data from this analysis as shown by the solid curves in Fig. |^. The fit yielded the following 
parameter values: < = 0.13 ± 0.03, < = 0.32 ± 0.14, 7 = 3.9 ± 0.9, ct = 0.11 ± 0.01, cl = 
0.12 ± 0.03, Cg = 0.014 ± 0.002; the errors include the statistical and systematic uncertainties 
combined in quadrature. 

The BEKW model also describes the dependence of shown by the solid curves 

in Fig. |13|, and the (3 (and Q'^) dependence of F.[^^'^^ shown in Fig. |ll|. The value 7 = 3.9 ±0.9 is 



consistent with the prediction of |]6^, 7 = 3, and somewhat higher than the value 7 = 2 given 
2^. 



m 



It is instructive to compare the (3 and dependences of the three components which build 
up the diffractive structure function F^^^^ in the BEKW model using the results from the fit. 
Figure |l4|(top) shows c-pF^ (dashed), clF^ (dashed-dotted), CgF^^ (dotted) and their sum 

XjpF^^^\xjp, (3,Q'^) at Xjp = Xo (solid curves) as a function of /3 for = 8,14,27,60 GeV^. 
Our data suggest that for /3 > 0.2 the colourless system couples predominantly to the quarks 
in the virtual photon. The region f3 > 0.8 is dominated by the contributions from longitudinal 
photons 0. The contribution from coupling of the colourless system to a qqg final state becomes 



important for P < 0.3. The last result is in contrast to the HI observation Jl^ that, using a 
DGLAP NLO fit, the large f3 region is dominated by the gluon contribution ^. 

Figure ^(bottom) shows the same quantities as a function of for j3 = 0.1,0.5,0.9. The 
gluon term, which dominates aX (3 = 0.1 rises with while the quark term, which is important 
at (3 = 0.5 shows no evolution with Q^. The contribution from longitudinal photons, which is 
higher twist and dominates at /3 = 0.9, decreases with Q"^. 

*In determining the diffractive cross section and the diffractive structure function the term 

[1 — ^dij/^^di// '] l^^s been neglected, see Eq. ^. If this term is kept for (3 > 0.8 and the BEKW 

fit is repeated with the assumption — <^t the changes in the fit parameters are small compared to their 
errors. 

^In two possible solutions were found from fits to the HI data: one where the gluon term is dominant 
at large f3 and one where it is not. The latter had a slightly larger value. 



In the BEKW model the x^p -dependence of the quark and gluon contributions for transverse 
photons is expected to be close to that given by the soft pomeron, tit ~ 2{oi^^°^^ ~ !)• However, 
perturbative admixtures in the diffractive final state are expected to have a somewhat stronger 
energy dependence, leading to an effective ut > 2(cr^*°'^* — 1). The dependence of the 
longitudinal contribution is driven by the square of the proton's gluon momentum density 
leading to 11^ > rij-. The fit results agree with these predictions but the errors are too large for 
a definitive statement. 

The same conclusion is reached when separate fits are performed for the regions (3 > 0.8 
and (3 < 0.8. Assuming n = rij- = the results are n{j3 > 0.8) = 0.46 ± 0.12 and 
n{p < 0.8) = 0.27 ± 0.03. It is important to note that already at = 8 GeV^ riTiQ"^ = 
8 GeV^) = 0.25± 0.04 which is substantially larger than the expectation for soft contributions, 
nsoft = 0.152tQ;o^8, indicating that the transverse and gluon components receive sizeable 
contributions from perturbative processes. 



In the BH model the and dependences of the diffractive structure function at small 
Xjp have been related to the x and dependences of the structure function F2 by assuming 
that in diffractive DIS a colourless cluster a is separated from the proton which interacts with 
the virtual photon. The probabihty, cr(x^,Q^), for finding such a cluster in the proton at 
small Xjp is expected to have an x^ dependence similar to the x dependence of the quark and 
gluon densities in the proton, g{x,Q^), qsea{x,Q^), provided x = x^. Since diffractive DIS is 
expected to predominantly produce configurations where the relative transverse momenta of at 
least one pair of partons are small, QCD evolution is suppressed in contrast to inclusive deep- 
inelastic scattering. These arguments have led to the prediction XjpF^^^\xjp, P^Q"^) oc F2{x = 
X p , Q"^) / logiQ (Q'^/Ql). Here, Ql ^ 0.55 GeV^ was taken from an analysis of the F2{x,Q'^) 
data from HERA |6T|. This relation predicts similar W dependences for the diffractive and total 
cross sections which is in agreement with the data presented above. It also predicts different 
dependences for x^F^^^' and F2. 



The behaviour of the two structure functions is compared in Fig. |T2] which shows 



Xp F^^^\xjp, PjQ"^) (solid points) for fixed values of Mx and W and F2{x = x^^Q"^) (open 
points). The F2{x = x^, Q"^) values were obtained from our published data [1^ by taking that 
measurement of ^2(0:, Q^) with x closest to x^ and transporting it to x = x^p. For ease of 
comparison, F2 has been multiplied by a constant factor of 0.06. The comparison shows that 
the two quantities have different evolution with Q"^. For Mx < 7.5 GeV, x^F^^^^ decreases 
with while the structure function F2{x,Q'^) gradually rises with Q^. In Fig. ^ x^F^^^^ 
is also compared with F2{x = x jp , Q^) / logiQ^Q"^ / Ql) (points marked as stars) as suggested 
by [^. Here, the F2 values were multiplied by a factor of 0.05 as obtained from a fit to the 
x^F2^^^ for Mx = 2 GeV. The evolution of the data at low Mx {Mx < 3 GeV) is well 
described by this model. At larger Mx values there is a tendency for the data to lie above 
the BH prediction. We note that in the BEKW model this is understood as resulting from the 
logarithmic growth of the {qqg) contribution with Q^. 



13 Summary and conclusion 

The DIS diffractive cross section da^fp_^xN/dMx, has been measured for Mx < 5.5 GeV, Mx < 
15 GeV, 60<W < 200 GeV and 7 < < 140 GeV^. For fixed the diffractive cross section 
rises rapidly with W. A fit of the W dependence by the form da'^Ip^xNi^x,W,Q'^)/dMx oc 



W"- " yielded a'^'^f = 0.507±0.034(stat)i|]:J^i (syst) which corresponds to a t-averaged pomeron 
trajectory of «^ = 1.127 ± .009 {stat) to fil (syst). The rise is faster than expected in Regge 
models using the intercept of the pomeron trajectory extracted from hadron-hadron scattering. 
The W dependence of the diffractive cross section, contrary to naive expectations, is the same 
as that of the total virtual photon proton cross section. The diffractive contribution to the 
total cross section for Mx < 15 GeV, Mjy < 5.5 GeV and ISA < W < 200 GeV amounts to 
13.2i°:^% at ^ 8 GeV^ decreasing to 3.7lgi% at = 60 GeV^. 

The analysis of the data in terms of the diffractive structure function F^^^\xjp, (3,Q'^) of the 
proton shows that XjpF^^^^ rises as 0. The data are consistent with the assumption that 

the diffractive structure function F^^^^ factorizes into a term depending only on x^ and a struc- 
ture function F^^'^^ which depends on ((3,Q'^). The rise of XjpF^^^^ with x^ can be described 

as XjpF2^^^ oc {l/x^Y with n = 0.253 ± 0. 017 {stat)to^ll{syst). The data are also consistent 
with models which break factorization. The rise of ^2^*''^'' reflects the rise of da'!^!p^xN / 

with W . For flxed Mx < 7.5 GeV and flxed H^, x^F^^^^ decreases slowly with Q'^ while for 
Mx = 11 GeV it is approximately constant. 

The data have been compared with several partonic models of diffraction. Good agreement with 
the data can be achieved. The models provide a flrst glimpse of how the different components 
may build up the diffractive structure function. The behaviour of XpF^^^\xjp, l3,Q'^) is 
different from that of the proton structure function F2{x,Q'^), taken at x — x^, which rises 
gradually with Q'^. It is in broad agreement with the BH conjecture that XjpF^^^\xjp, Q^) oc 
F2{x = x^,g2)/iog^p(Q2/g2) where Ql = 0.55 GeV^. Using the BEKW model at medium 
j3 the main contribution comes from transverse photons coupling to a qq system. The region 
P < 0.2 is dominated by qqg contributions. Longitudinal photons coupling to a qq system 
account for most of the data at > 0.8. The transverse photon qq contribution, which is 
dominant, is of leading twist and has no substantial evolution with Q^. 

The leading twist behaviour and the strong rise of XjpF^^^^ as a;^ — > suggest a partonic 
process as a major production mechanism for diffractive scattering in DIS. 
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Table 1: Cross section for diffractive scattering via 7*p XN, where is the proton or 
dissociated nucleonic system with mass < 5.5 GeV function of Mx,Q^ and W. The 
statistical and systematic errors are given. The overall normalization uncertainty of ±2% is 
not included. 



Mx 




H' 


da^VJ ,^^/dMx 


± Stat 


± SA'St 


(Ge\-) 


(Ge\-^') 


(GcV) 




(iilj/C;c\') 




2.0 
2.0 
2.0 
2.0 
2.0 


8.0 
8.0 
8.0 
8.0 
8.0 


66.7 
81.8 
99.8 
122.1 
148.6 


138.6 
153.8 
180.0 
192.3 
221.9 


± 9.5 
± 10.1 
± 11.7 
± 13.4 
± 13.2 


+ 11. S 
— 11.6 
+ 11.8 
— 12.2 
+ 13.6 
-11.7 
+ 19.6 
-13.3 
+20.2 
-16.5 


2.0 


8.0 


181.5 


219.9 


± 13.9 


+22.2 
— 17.5 


5.0 
5.0 
5.0 
5.0 


8.0 
8.0 
8.0 
8.0 


66.9 
81.6 
99.9 
121.9 


120.3 
127.9 
139.8 
162.8 


± 7.9 
± 7.0 
± 7.5 
± 8.2 


+ 11. b 

— 27. 1 
+ 13.1 

-21.2 
+ 12.1 
-17.7 
+8.0 

— 14.0 


5.0 
5.0 


8.0 
8.0 


148.9 
181.7 


175.4 
198.4 


± 8.8 
± 10.0 


+20.6 

— 12.4 
+11.0 

— 13.6 


11.0 

11.0 
11.0 


8.0 

8.0 
8.0 


121.7 

149.1 
181.6 


68.9 
74.5 
78.4 


± 4.2 

± 4.4 
± 4.8 


+ 5.9 

— 22.9 
+ 7.0 

— IG.l 
+ 2.7 


2.0 


14.0 


67.2 


48.0 


± 3.6 


+3.V 
— 6.1 


2.0 
2.0 
2.0 
2.0 


14.0 
14.0 
14.0 
14.0 


81.5 
100.0 

121.6 
148.5 


48.9 
58.1 
60.0 
55.2 


± 4.1 
± 4.8 

± 5.5 
± 5.7 


+9.2 
—4.2 
+7.2 
—4.0 

+ 10.6 

— 4.4 
+ 15.6 

— 5.0 


2.0 


14.0 


181.8 


76.7 


± 5.3 


+6.9 


5.0 
5.0 

5.0 


14.0 
14.0 
14.0 


67.0 
81.7 

99.9 


52.3 
61.4 
73.5 


± 3.7 
± 3.2 
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— 12.7 
+ 5.5 
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+ 3.8 
— 9 . 5 
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+9.1 
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83.8 
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± 2.4 


+1:7 
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27.0 
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13.4 


± 2.0 
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27.0 
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± 2.3 
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27.0 


121.4 
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+ 3.4 


2.0 


27.0 


148.9 


20.5 


± 3.2 


+ 2^8 
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27.0 


81.7 
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26.4 
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149.2 
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— 31 
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27.0 
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33.4 


± 2.7 


+4:1 
-2.7 
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27.0 


148.9 


19.6 


± 1.6 


+ 1.6 
-4.4 


11.0 


27.0 


181.5 


25.8 


± 2.0 


+0.6 
-4.1 


2.0 


60.0 


81.2 


0.8 


± 0.3 


+0.6 
-0.3 


2.0 
2.0 
2.0 


60.0 

60.0 
60.0 


101.1 

122.5 
148.8 


1.9 
1.4 
2.1 


± 0.5 

± 0.4 
± 0.5 
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-0.6 
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-0.2 
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± 1.0 
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122.5 


4.4 


± 0.6 
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3.9 
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+ 1.5 
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5.0 


60.0 
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6.1 


± 0.8 


+ 1.2 
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11.0 


60.0 
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± 0.8 


+ 1.U 
-1.7 


11.0 
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181.6 


7.8 


± 0.8 


+0.3 
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Table 2: The diffractive structure function multiplied by x^, XjpF^^^\xjp, (3, Q^), for diffractive 
scattering via XN, where is a nucleonic system with mass Mjy < 5.5 GeV as a function 

of Xjp, j3 and Q^. The statistical and systematic errors are given. The overall normalization 
uncertainty of ±2% is not inchidcd. 



(GcV^) 



''IP ^2 



± stat ± syst 



-0.0025 
+0.0025 
-0.0026 
+0.0029 
-0.0025 
+0.0042 
-0.0028 
+0.0043 
-0.0035 
+0.0047 
-0.0037 
+0.0027 
-0.0064 
+0.0031 
-0.0050 
+0.0028 
-0.0042 
+0.0019 
-0.0033 
+0.0049 
-0.0029 
+0.0026 
-0.0032 
+0.0025 
-0.0096 
+0.0029 
-0.0067 
+0.0011 
-0.0054 



0.00269 
0.00179 
0.00120 

0.00081 
0.00054 
0.00036 
0.00735 
0.00495 
0.00330 
0.00222 
0.00149 
0.00100 
0.00871 
0.00580 
0.00391 



0.667 
0.667 
0.667 
0.667 
0.667 
0.667 
0.242 
0.242 
0.242 
0.242 
0.242 
0.242 
0.062 
0.062 
0.062 



8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 



0.0297 
0.0329 
0.0385 
0.0411 
0.0475 
0.0471 
0.0284 
0.0301 
0.0329 
0.0383 
0.0413 
0.0467 
0.0288 
0.0312 
0.0328 



± 0.0020 
± 0.0022 
± 0.0025 
± 0.0029 
± 0.0028 
± 0.0030 
± 0.0019 
± 0.0016 
± 0.0018 
± 0.0019 
± 0.0021 
± 0.0023 
± 0.0018 
± 0.0018 
± 0.0020 



+11.01)21 
-0.0034 
+ 0.0052 
-0.0023 
+0.0041 
-0.0022 
+0.0060 
-0.0025 
+0.0087 
-0.0028 
+0.0039 
-0.0040 
+0.0025 
-0.0062 
+0.0027 
-0.0058 
+0.0018 
-0.0046 
+0.0023 
-0.0036 
+0.0044 
-0.0029 
+0.0059 
-0.0016 
+0.0028 
-0.0074 
■ 0-()()28 
-0.0(148 



0.00398 
0.00270 
0.00180 
0.00122 
0.00082 
0.00054 
0.00866 
0.00583 
0.00390 
0.00263 
0.00176 
0.00118 
0.00607 
0.00109 



0.778 
0.778 
0.778 
0.778 
0.778 
0.778 
0.359 
0.359 
0.359 
0.359 
0.359 
0.359 
O.lO'l 
0.101 



14.0 
14.0 
14.0 
14.0 
14.0 
14.0 
14.0 
14.0 
14.0 
14.0 
14.0 
14.0 
14.0 
1 1.0 



0.0270 

0.0275 
0.0327 
0.0337 
0.0310 
0.0431 
0.0255 
0.0300 
0.0358 
0.0375 
0.0408 
0.0408 
0.0306 
O.O:-!:-!! 



± 0.0020 
± 0.0023 
± 0.0027 
± 0.0031 
± 0.0032 
± 0.0030 
± 0.0018 
± 0.0016 
± 0.0018 
± 0.0018 
± 0.0020 
± 0.0021 
± 0.0017 
± O.OOIS 



+(1.0038 
-0.0033 
+0.0052 
-0.0032 
+0.0048 
-0.0024 
+0.0063 
-0.0031 
+0.0053 
-0.0050 
+0.0031 
-0.0063 
+0.0040 
-0.0068 
+0.0025 
-0.0045 
+0.0041 
-0.0037 
+0.0020 
-0.0039 
+0.0051 
-0.0034 
+0.0026 
-0.0071 
+0.0010 
-0.0067 



0.00682 
0.00457 
0.00313 
0.00210 
0.00140 
0.00093 
0.00776 
0.00524 
0.00350 
0.00233 
0.00158 
0.00667 
0.00449 



0.871 
0.871 
0.871 
0.871 
0.871 
0.871 
0.519 
0.519 
0.519 
0.519 
0.519 
0.182 
0.182 



27.0 
27.0 
27.0 
27.0 
27.0 
27.0 
27.0 
27.0 
27.0 
27.0 
27.0 
27.0 
27.0 



0.0172 
0.0251 
0.0233 
0.0299 
0.0382 
0.0453 
0.0265 
0.0296 
0.0331 
0.0411 
0.0418 
0.0318 
0.0417 



± 0.0027 
± 0.0037 
± 0.0043 
± 0.0049 
± 0.0059 
± 0.0066 
± 0.0026 
± 0.0025 
± 0.0029 
± 0.0032 
± 0.0034 
± 0.0026 
± 0.0032 



+0.0052 
-0.0026 
+0.0017 
-0.0053 
+0.0033 
-0.0021 
+0.0106 
-0.0027 
+0.0035 
-0.0230 
+0.0031 
-0.0079 
+0.0045 
-0.0049 
+0.0068 
-0.0026 
+0.0055 
-0.0028 
+0.0043 
-0.0074 
+0.0014 
-0.0057 



0.00961 
0.00622 

0.00425 
0.00288 
0.00197 
0.00856 
0.00564 
0.00381 
0.00255 
0.00817 
0.00548 



0.938 
0.938 
0.938 
0.938 
0.938 
0.706 
0.706 
0.706 
0.706 
0.331 
0.331 



60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 



0.0073 
0.0167 

0.0119 
0.0178 
0.0386 
0.0191 
0.0202 
0.0179 
0.0276 
0.0258 
0.0345 



± 0.0023 
± 0.0045 
± 0.0032 
± 0.0046 
± 0.0089 
± 0.0030 
± 0.0029 
± 0.0032 
± 0.0036 
± 0.0034 
± 0.0036 
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Figure 1: Reaction 7*p X + anything, where X is the system observed in the detector. Top: 
Distributions of Mx, the corrected mass of the system X. The distributions are not corrected 
for acceptance effects. The shaded histograms show the distributions of events with rjmax < 1-5. 
Bottom: Same distributions as above presented in terms of lnM|-. The straight hues give the 
nondiffractive contributions as obtained from the fits. The upper curves show the fit resuhs for 
the sum of the diffractive and nondiffractive contributions. 




Figure 2: Reaction 7*^9 ^ X + anything, where X is the system observed in the detector. Top: 
Distributions in lnM|^ as predicted by ARIADNE for the nondiffractive contribution at the 
generator level (solid histograms) and detector level (dashed histograms) for the W intervals 
and values indicated. The straight lines show the results of the fits to the distributions at the 
detector level. Bottom: Distributions in InM^ for the sum of the diffractive and nondiffractive 
contributions as predicted at the detector level by RAPGAP plus ARIADNE (points with error 
bars) and for the nondiffractive contribution alone (dashed histograms). The straight lines show 
the results for the nondiffractive contribution obtained from fitting the sum of the diffractive 
and nondiffractive contributions with fe^f,^. The upper curves show the fit results for the sum 
of the diffractive and nondiffractive contributions. 



W = 134-164 GeV 
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Figure 3: Reaction 7*p X +anything, where X is the system observed in the detector. Shown 
are distributions of lnM|- at the detector leveL The dotted histograms show the diffractive 
contributions from 7*p — > XA^*^**"^ as predicted by EPSOFT. The sohd histograms show the 
sum of the diffractive contributions from 7*p — > XN^^^^°'^ and 7*p — > Xp (the latter as predicted 
by RAPGAP). The dashed histograms show the nondiffractive contributions as predicted by 
ARIADNE. The points with error bars show the sum of the diffractive and the nondiffractive 
contributions. The dashed straight hues show the fits performed to the nondiffractive contribu- 
tions alone. The straight lines show the results for the nondiffractive contribution from fitting 
the sum of the diffractive and nondiffractive contributions with . 




Figure 4: Reaction 7*p — > X + anything, where X is the system observed in the detector. 
Distributions in lnM|^ for data. The straight hues give the nondiffractive contributions and 
the upper curves the sum of the diffractive and nondiffractive contributions as obtained from 
the fits. The sohd histograms show the distributions for events with a PRT signaL 
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Figure 5: The differential cross sections da^Ip^xN/dMx, < 5.5 GeV, as a function of 
W at average values of Mx = 2, 5, 11 GeV, = 8, 13, 27, 60 GeV^. The inner error 
bars show the statistical errors and the full bars the statistical and systematic errors added 
in quadrature. The overall normalization uncertainty of 2% is not included. The solid curves 
show the result from fitting the diffractive cross section for each {W, Q^) bin separately using 
the form da^fp^xN/^^x oc (ly^)"*'''' where a'^*-^-^ and the normalization constants were treated 
as free parameters. The dashed curves show the result from the fit where a'^^-^^ was assumed to 
be the same for all [W, Q"^) bins. 
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Figure 6: The diffractive differential cross section for 7*p — > XN, Mjv < 5.5 GeV, multiplied 
by Q^, Q'^da'^'-ff {YP XN)/dMx, as a function of Q"^ for the values of Mx and W indicated. 
The inner error bars show the statistical errors and the full bars the statistical and systematic 
errors added in quadrature. 
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Figure 7: The parameter cF^ obtained from the fits versus for Mx = 2 and 5 GeV. The 
inner error bars show the statistical errors and the full bars the statistical and systematic errors 
added in quadrature. The result for a soft pomeron is indicated by the heavy horizontal band. 
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Figure 8: The ratio of the diffractive cross section, integrated over the Mx intervals indicated, 



diff 



dMx(j'^*p^xN^ fo^^ < 5.5 GeV, to the total cross section for virtual photon 
proton scattering, rfj^ = a'^'-^-^/cr^^p, as a function of W for the Mx intervals and values 
indicated, cj^o/ was taken from our F2 measurements using the 1994 data [jl9|. The inner error 
bars show the statistical errors and the full bars the statistical and systematic errors added in 
quadrature. 
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Figure 9: The diffractive structure function of the proton multiphed by x^p, XjpF^^^\ as a 
function of Xjp from this analysis (solid points). The inner error bars show the statistical errors 
and the full bars the statistical and systematic errors added in quadrature. The curves show the 
results from the models of Nikolaev and Zhakarov (NZ), Bialas, Peschanski and Royon (BPR) 
and Bartels, Ellis, Kowalski and Wiisthoff (BEKW). 
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Figure 10: The diffractive structure function of the proton for 7*p — > XN, Mj^ < 5.5 GeV, 
multipHed by x^, XjpF^^^\xjp, P,Q'^), from this analysis (sohd points) compared with the 
results from our previous LPS measurement obtained with an identified proton for 7*p — >■ Xp 
(stars) and from a subsample of the HI data (open points) for 7*p XN, Mjy < 1.6 GcV. For 
ease of comparison the results from this analysis are scaled to {f3,Q^) values used by HI. The 
data points from this experiment shown for = 7.5 and 9 GeV^ are those obtained at = 8 
GeV^ shifted to = 7.5 and 9 GeV^. The LPS data are given for slightly different P and 
values. 
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Figure 11: The structure function Fp^'{(3,Q^) for YP XN,Mn < 5.5 GeV, for the 
values indicated, as a function of (3 as extracted from a fit to the measured x^F^^^^^^ values, 
see text. The error bars show the statistical and systematic errors added in quadrature. The 
curves show the fit results obtained with the BEKW model discussed in the text. 
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Figure 12: The solid points show the diffractive structure function for 7*p —>■ XN, Mjy < 
5.5 GeV, multiphed by Xjp, XjpF^^^\xjp, (3,Q'^), from this analysis, as a function of for 
the {Mx,W) values indicated. The inner error bars show the statistical errors and the full 
bars the statistical and systematic errors added in quadrature. The open points show the 
structure function F2{x = Xjp,Q'^) multiplied by 0.06 obtained with the 1994 data [|l^ as a 
function of for the {MxyQ"^) values indicated. The points marked by stars show F2{x = 
Xjp, Q^) / logioiQ"^ / Qq) scaled by a factor 0.05. Note that Xjp is known when Mx, W and are 
given. 
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Figure 13: The diffractive structure function of the proton for 7*p — > XN, Mj^ < 5.5 GeV, 
multiphed by x^, XjpF^^^\xjp, (3,Q'^), as a function of from this analysis. The inner error 
bars show the statistical errors and the full bars the statistical and systematic errors added in 
quadrature. The curves show the fit results obtained with the BEKW model. 
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Figure 14: Top: The three components {qq)T, i^Qg) and {qq)L of the BEKW model building 
up the diffractive structure function of the proton and their sum F<^^'^\j3 , Q"^) as a function 
of p for Q2 = 8, 14, 27 and 60 GeV^, as obtained from a fit of the BEKW model to the data. 
Bottom: the same quantities as a function of for /3 = 0.1, 0.5 and 0.9. 



